LARGE, GLOBAL SOLUTIONS TO THE NAVIER-STOKES EQUATIONS, 
SLOWLY VARYING IN ONE DIRECTION 



JEAN- YVES CHEMIN AND ISABELLE GALLAGHER 



Abstract. In [3] and [4] classes of initial data to the three dimensional, incompressible 
Navier-Stokes equations were presented, generating a global smooth solution although the 
norm of the initial data may be chosen arbitrarily large. The aim of this article is to provide 
new examples of arbitrarily large initial data giving rise to global solutions, in the whole space. 
Contrary to the previous examples, the initial data has no particular oscillatory properties, 
but varies slowly in one direction. The proof uses the special structure of the nonlinear term 
of the equation. 



1. Introduction 

The purpose of this paper is to use the special structure of the tridimensional Navier-Stokes 
equations to prove the global existence of smooth solutions for a class of (large) initial data 
which are slowly varying in one direction. Before entering further in the details, let us recall 
briefly some classical facts on the global wellposedness of the incompressible Navier-Stokes 
equations in the whole space M'^. The equation itself writes 

{dtu + u ■ Vu — Au = —Vp 
div M = 
uit=o = uo 

where u = = {u^,u^) is a time dependent vector field on M^. The divergence free 

condition determines p through the relation 

-Ap= djdkiu^u''). 
This relation allows to put the system {NS) under the more general form 



(GNS) I - = Q{u, u) 
1 u\t=o = Uo 

where Q{v,w) Qj^k{D){v^w^) and Qj^k{D) are smooth homogeneous Fourier multi- 

i<i,fc<3 

pliers of order 1. 

Moreover, this system has the following scaling invariance: if {u,p) is a solution on the time 
interval [0,r), then {ux,px) defined by 

u\{t, x) Xu{}?t, Ax) and P\{ti x) }?p[}?t, Ax) 



Key words and phrases. Navier-Stokes equations, global wellposedness. 
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is a solution on the time interval [0, A~^T). Of course, any smallness condition on the initial 
data that ensures global solutions, must be invariant under the above scaling transformation. 
The search of the "best" smallness condition is a long story initiated in the seminal paper of 
J. Leray (see |12j). continuated in particular by H. Fujita and T. Kato in [6], Y. Giga and 
T. Miyakawa in [8], and M. Cannone, Y. Meyer and F. Planchon in [T]. This leads to the 
following theorem proved by H. Koch and D. Tataru in [llj . In the statement of the theorem, 
P{x, R) stands for the parabolic set [0, R'^] x B{x, R) where B{x, R) is the ball centered at x, 
of radius R. 

Theorem 1 ([llj). If the initial data uq is such that 

def 1 f 

(1.1) IkollljMO-i = supt||e*^no||i^ + sup / |(e*^uo)(t, y)pd?/ 

t>0 xeRS -K Jp{x,R) 

R>0 

is small enough, then there exists a global smooth solution to {GNS). 



A typical example of application of this theorem is the initial data 

(1.2) Uo(x) =^cos^^^(92(/>(xi,X2),-9i0(2;i,X2),O) with 

as soon as e is small enough (see for example [3] for a proof). The above theorem is probably 

the end point for the following reason, as observed for instance in [4|. If 2? is a Banach 

space continuously included in the space <S' of tempered distributions on M^, such that, for 

any (A, a) G M+ x R^, ||/(A(.-a))||ij = A'^ then \\-\\b < Csupt^ ||e*^tio||Loc. The second 

t>o 

condition entering in the definition of the BMO^^ norm given in (II. Ih merely translates the 
fact that the first Picard iterate should be locally square integrable in space and time. 

Those results of global existence under a smallness condition do not use the special structure 
of the incompressible Navier-Stokes system and are valid for the larger class of systems of the 
type iGNS). The purpose of this paper is to provide a class of examples of large initial data 
which give rise to global smooth solutions for the system {NS) itself, and not for the larger 
class (GNS). In all that follows, an initial data uq will be said "large" if 

/I o\ II II def , 1 II tA II 

(l.oj II^oIIr-1 =supt2||e uq\\l^ 

oo,oo 

is not small. 

Such initial data, in the spirit of the example provided by (|1.2p . are exhibited in ^ (see 
also [3] for the periodic case). In particular, the following theorem is proved in [4j. 

Theorem 2 ([4j). Let G 5(M^) be a given function, and consider two real numbers e and a 
in ]0, 1[. Define 

^ {x) = cos[—)<t)[xi,—,x^y 

Then for e small enough, the smooth, divergence free vector field 

Ul{x) = (5299^(2;), -9193^ (x),0) 

satisfies lirn ||tioll^-i = 00, and generates a unique global solution to (NS). 



LARGE, GLOBAL SOLUTIONS TO THE NAVIER-STOKES EQUATIONS 



3 



The proof of that theorem shows that the solution remains close to the solution of the free 
equation e*^uo,e. It is important to notice that the proof uses in a crucial way the algebraic 
structure of the non linear term u ■ Vti, but uses neither the energy estimate, nor the fact 
that the two dimensional, incompressible Navier-Stokes is globally wellposed. Let us give 
some other results on large initial data giving rise to global solutions: In [13], the initial 
data is chosen so as to transform the equation into a rotating fluid equation. In [lOj . |14j 
the equations are posed in a thin domain (in all those cases the global wellposedness of the 
two dimensional equation is a crucial ingredient in the proof). In p[j and the case of an 
initial data close to bidimensional vector field is studied, in the periodic case. Finally in [5] 
(Remark 7), an arbitrary large initial data is constructed in the periodic case, generating a 
global solution (which is in fact a solution to the heat equation, as the special dependence on 
the space variables implies that the nonlinear term cancels). 

The class of examples we exhibit here is quite different. They are close to a two dimensional 
flow in the sense that they are slowly varying in the one direction (the vertical one). More 
precisely the aim of this paper is the proof of the following theorem. 

Theorem 3. Let Vq = {vq,Vq,0) be a horizontal, smooth divergence free vector field on M.^ 
(i.e. Vq is in L^(IR^) as well as all its derivatives), belonging, as well as all its derivatives, 
to L^(Ma,3; H~^{M.'^)); let wq be a smooth divergence free vector field on M^. Then if e is small 
enough, the initial data 

Uo(x) = {vq + ew'',Wo){xh,ex3) 
generates a unique, global solution u'^ of (NS). 

Remarks 

• A typical example of vector fields Vq satisfying the hypothesis is Vq = {—d24>,di(t),0) 
where (f) is a function of the Schwarz class 5(M^). 

• This class of examples of initial data corresponds to a "well prepared" case. The "ill 
prepared" case would correspond to the case when the horizontal divergence of the 
initial data is of size e"^ with a less than 1, and the vertical component of the initial 
data is of size e°^~^. This case is certainly very interesting to understand, but that 
goes probably far beyond the methods used in this paper. 

• We have to check that the initial data may be large. This is ensured by the following 
proposition. 

Proposition 1.1. Let {f,g) be in 5(M^). Let us define h^{xh,X3) f{xh)gisx3). We have, 
if e is small enough. 




Proof. By the definition of || • ^j^a^ given by (|1.3p . we have to bound from below the 

quantity ||e*^^^||Loo(]83-). Let us write that 

(e*^/i^)(t,x) = (e''"^f){t,Xh){e'^lg){eH,ex3). 
Let us consider a positive time Iq such that 

4l|e*''^'^/llLoc(M^)>^ll/b-^(M^). 
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Then we have 

4l|e*«^/i^|li^(M3) = 4l|e*''^VIlLoo(K2)||(e*«^3<7)(e'to,e-)llLoo(M) 

^ ^ll/llB-^{M^)l|e''*°'''5llL-(R)- 

As hm ^od^^g _ ^ L°°(M), the proposition is proved. □ 

Structure of the paper: The proof of Theorem [3] is achieved in the next section, assuming 
two crucial lemmas. The proof of those lemmas is postponed to Sections [3] and H] respectively. 

Notation: If A and B are two real numbers, we shall write A < 5 if there is a universal 
constant C, which does not depend on varying parameters of the problem, such that A < CB. 
If ^ < B and B <A, then we shall write A ~ 5. 

If fo is a vector field, then we shall denote by a constant depending only on norms of vq. 
Similarly we shall use the notation C^o,«;o if the constant depends on norms of two vector 
fields vq and wq^ etc. 

A function space with a subscript "/i" (for "horizontal") will denote a space defined on M^, 

while the subscript "w" (for "vertical" ) will denote a space defined on M. For instance ^ 

LP{m.^), Ll = L''(M), and similarly for Sobolev spaces or for mixed spaces such as L^Vj^ 
or LlHl. 



2. Proof of the theorem 

The proof of Theorem [3] consists in constructing an approximate solution to {NS) as a per- 
turbation to the 2D Navier-Stokes system. Following the idea that we are close to the two 
dimensional, periodic incompressible Navier-Stokes system, let us define as the solution of 
the following system, where ys G M is a parameter: 

{dtv^ + y!" ■ Vhv!" - Ahv!" = -VhPQ in M+ x 
divhn'' = 
l^S=o = ^o(->y3)- 

This system is globally wellposed for any € M, and the solution is smooth in (two dimen- 
sional) space, and in time. Let us consider the solution uf of the linear equation 

r dtwf + ■ Vhvf - Ahvf - e'^dluf = -{V^PvE^d^ip^) in M+ x 

(r|) I diYvf = 

and let us define the approximate solution 

(2.1) vl^{t,x) = ((^;^0) + e(u;^•^e-V'^))(^,Xh,ex3) and 

Plppit^x) = {p^ + ep^)it,Xh,ex3). 

Finally let us consider the unique smooth solution of (NS) associated with the initial 
data Uq on its maximal time interval of existence [0, T^). The proof of Theorem [3] consists in 
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proving global in time estimates on v^pp, in order to prove that i?' 
globally in time; this ensures the global regularity for {NS). 



def 



— Vapp remains small, 



More precisely, the proof of Theorem [3] relies on the following two lemmas, whose proofs are 
postponed to Sections [3] and H] respectively. 



Lemma 2.1. The vector field v^pp defined in i2.1\) satisfies the following estimate: 



\v, 



app 1 1 L2 ( 



DO ,wo ■ 



Lemma 2.2. The vector field W 



def 



u 



•^app 



satisfies the equation 



dtR" + ■ VR' - AR' + ■ VR^ + R' ■ Vvl 



(En 



app 

div R^ 







^app 



R\t=0 = 



with \\F^\ 



'L2( 



1 



<a 



vo,wo 



1 

£3 . 



Let us postpone the proof of those lemmas and conclude the proof of Theorem [3l We denote, 
for any positive A, 

^ onrl ??? 'i^* o^T. / _ \ f 



ydt) =^ \K„jt, 



and Rl{t)=^ exp(^-X J^Ve{t')dt'^ R%t). 



Lemma l2 . 1 1 implies that Iq 
1 d 



def 



Vs{t)dt is finite. By an H2 energy estimate in M^, we get 



2dt' 
+ 



■ 1 + l|Vi?i 



\ 1 < -2AK 



{RUt)-VRUt)\RUt))^r 



{Ri{t)-Vvlpp{t)\Rm 



+ 



{v'{t)-Rm\Rx{t)) 



The estimate (i) of Lemma 1.1 of [2j claims in particular that, for any s g] — (i/2,d/2[, for 
any divergence free vector field a in d space dimensions and any function 6, we have 



(2.2) 



(a.V6|%. <C||Va||^._J|6||^.||V6||^.. 



xWVRl 



Applying with d = 2> and s = 1/2, this implies that 
(2.3) |(^i(t)-Vi?i(t)|i?i(t))^i|<||i, 
In order to estimate the other non linear terms, let us establish the following lemma. 
Lemma 2.3. Let a and b he two vector fields. We have 



2 

. 1 . 
_f/2 



(a • Vb\b) 



+ 



{b-Va\b)^. 



\L^ + \\Va\\L^^Ll)) l|fo||^i||V6||^i. 



Proof. By definition of the H2 scalar product, we have 

(a-V6|6)^i < \\a-Vbh2\\Vbh2 
< l|a||L-||V6||i2. 
The interpolation inequality between Sobolev norm gives 



{a-Vb\b)^, < llallz.^ 11611^1 II V6|| 
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Now let US estimate (6 ■ Va|6) . i . Again we use that 



Then let us write that 



\b-Va\\\.2 = / \b{xh,xz)\Ia{xh,x:>,)\'^dxhdxz. 



Gagliardo-Nirenberg's inequality in the horizontal variable implies that 



VX3 GM, |6(X;„X3)P < ||6(-,X3)|| . i||V/,6(-,X3)|| . 1. 

Let us use the Cauchy-Schwarz inequality; this gives 

\\b-Va\\l2 < [ \\b{-,X3)\\ i\\VhK-,X3)\\ i\\Va{;X3)\\l2dx3 



< i|V«lli-(L2) / ||6(-,X3)|| i||V;,6(-,X3)|| idX3 



When s is positive, we have, thanks to Fourier-Plancherel in the vertical variable, 



/ / \^h\^'\J'hb{^h,X3)\^d^hdxs 
/ / m^-'mh,C3)fdChd^3 

JR JR2 



< 



< \\b\\% 



This concludes the proof of Lemma 12. 31 

Conclusion of the proof of Theorem [3l We infer from the above lemma that 
(i?i(t)-V<pp(OIi?i(t))^i + {v%,{t)-Rl{t)\Ri{t))^r 



<\\WRm\\\r+cvM\mt)\\\^^ 



Together with (|2.3p . this gives 



Choosing A such that C — 2A is negative, we infer that 

l^iWII^ + (1 - Ce^'«)||Vfli(t)i|^ < . 
Since i?^(0) = 0, we get, as long as . i is less or equal to l/2Ce~'*'^°, that 



\Rm\\\^^ + \j^ \\^RW)\\\.dt' <c. 
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We therefore obtain that W goes to zero in L°°(R+; ifa) nL2(M+; H^). That imphes that li^ 
remains close for ah times to f^pp, which in particular implies Theorem [3l □ 



3. Estimates on the approximate solution 

In this section we shall prove Lemma 12.11 stated in the previous section. The proof of the 
lemma is achieved by obtaining estimates on stated in the next lemma, as well as on nf 
(see Lemma 13.21 below) . 

Lemma 3.1. Let he a solution of the system {NS2Ds). Then, for any s greater than — 1 
and any q G N'^, we have, for any ys in M and for any positive t, 

-,2/3)11^. + r WVhv''{t',;y3)\\lJt' < a„(y3), 

h Jq h 

where C^o(') belongs to n and its norm is controled by a constant C^q. 

Proof. For s = and a = 0, the estimate is simply the classical L^-energy estimate with 
as a parameter: writing v = {v^, 0) we have 

(3.1) \\v{t,;y3)\\l2+2 f \\yhv{t',;y3)\\l2dt' = \\vo{;y3)\\l,. 

h Jq h h 

In the case when a = 0, the estimate (i) in Lemma 1.1 of [2] gives, for any s greater than — 1, 

~\\v{tr,y3)\\Hs^ + \\'^hv{tr,y3)\\%s^<c\\Vh^^^^ 

We infer that 

^Mt,;y3)\\ls + \\Vhv{t,;y3)\\ls < C\\Vhvit,;y3)\\l2\\vit,;y3)\\l^. 

dZ h h h h 

Gronwall's lemma ensures that 

Mt,;y3)\\%s+ I \\^hv{t\-^y3)\\\sdt' <\\vQ{-,y3)\\%e^v(c I \\Vhv{t' ,-,y3)\\l2dt'). 

h Jq h h \ Jq h / 

The energy estimate (j3.ip implies that 

\W,-,y3)f^^^ [ \\Vhvit\-,y3)\\lsdt' <\\M-,y3)\\%exp(c\\vQ\ 



This proves the lemma in the case when a = 0. Let us now turn to the general case, by 
induction on the length of a. It is clear that in the proof, we can restrict ourselves to the case 
when s s] — 1, 1[. 

Let us assume that, for some A; € N, 

(3.2) VsG]-l,l[, V (\\d''v{t,.,y3)\\%+ t\\d''Vhv{t',;y3)\\lsdt') <Ck,.M, 
with Ck,voi-) G L^nL'^iR). 
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Thanks to the Leibnitz formula we have, for |q;| < + 1, 

Performing a Hf^ energy estimate in the horizontal variable and using the estimate (j2.2p in 
the case when d = 2 gives 



2dV 



< C\\Vhv{t, •, ya) h'f Wd'^vit, yg)!!^. II V/,a"i;(t, -,^3)11 



0<a 



Hi 



To estimate the last term, we shall treat differently the case |/3| = and |/3| ^ 0. In the first 
case, we notice first that when s = 0, laws of product for Sobolev spaces in give 



<\\d>Htr,y3)\\Ll\\^hd>Ht,-,y3)h2\\VhvHtr,y3)\\Ll- 

h. h ti 

If s > 0, then again laws of product for Sobolev spaces in give, for s g]0, 1[, 



L2 



5V(t,-,y3) • V,2;'^(t,-,y3) a"i;(t,-,y3) l ^ < ||5V||,^. II V^i;''||^2 1| V^SVll^^ , 

h 

whereas if — 1 < s < 0, 

(5V(t,-,y3) • V;,/(t,-,y3)|5°i;(i,-,y3)) . < || V,,9V|Ih= II V;,i;^||i2 ||9V|Ih« • 

h 

So in any case we have 

(5V(t,-,y3) • V,/(t,-,y3)|9"i;(t,-,y3)) . < ^llV.aVll^. + C||V,i;"||2, ||9V||^.. 

Now let us consider the case when |/3j ^ 0. As the horizontal divergence of v is identically 0, 
we have 

(d'^-fv^-Wh 5V|5"2^(t,-,y3)) . < ||9'^-V(f,-,y3)09V(i,-,y3)||^.||V;,9°i;(i,-,y3)||^.. 
Laws of product for Sobolev spaces in give, for s g] — 1, 1[, 

l|5"-V®aV||^. < c||a°-V||^,,||9^V/,/||^._,,, 

where s' is chosen so that s < s' < 1. 
Finally we deduce that 



d 



^l|5Mi,-,y3)||^. + ||V,9°i;(t,-,y3)||^. < CUV.^zlli^ ||a"^;||^^ 



2^ pa„,ll2 



+ Ca V] ||5"-^(t,-,y3)||^.J9^V;,i;(t,.,y3)|U_jV;,9"i;(t,-,y3) 

I3<a 
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Gronwall's lemma together with the induction hypothesis (I3.2|) impUes that 



<( E I|9%(^y3)||^. +Cfc,,„(y3))exp(Cfc /"V/.2!(t',-,y3)|li2di'' 

^|Q|=fc + l ^ J Jo ^ 

The energy estimate (13. ip allows to conclude the proof of Lemma 13.11 □ 
From this lemma, we deduce the following corollary. 

Corollary 3.1. Let -i/^ he a solution of the system {NS2D3). Then, for any non negative a, 
we have 



Proof. To start with, let us assume a > 0. Lemma 13.11 applied with s = a — 1 implies that 
(3.3) > , Va G N , WA\l^(w+;LIhi) ^ ^^o- 

Then, for any non negative a, we have 



id3dXt,;y'3Wvit,;y'3)H.dy'^ 

-00 



< 2\\dsd'^v'it,-)h2H4d>\t,.)h.H. 
By the Cauchy-Schwarz inequality, we have 

- a h - 

From (13.31). we infer 



(3.4) Vct > , Va G N , ||5VllL2(K+;Looij;^) < C^o- 

Now, by interpolation, it is enough to prove the first inequality with a = 0. The sys- 
tem {NS2D3) can be written 

I ^'^-^'^={ with /t^ V Q,,{D){v^v'). 

where Qj^k are homogenenous smooth Fourier multipliers of order 1. By Sobolev embeddings 
in M^, we get, for any 7/3 in M, 

M-,y3)\\L^R+ xR^) < ll2^o(-'y3)llij^i + ll/(-.y3)llLi(M+;if-i) 

< \\vo{;y3)\\H-^+C\\vi;y3)f ^ .1 



As sup,„ ||ff-,?/3)l| 1 < llf II 1 , we infer from J3l 

Il!llr2m+ .Tt^3^ < |knll?.£r-i + C^nlll^lP i < Q 



— L2fM+ X ^ — r2 tr-i ~r "-^I'O ii , .1 



The corollary is proved. □ 
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Finally we have the following estimate on uf. 

Lemma 3.2. Let uf be a solution of the system (T^). Then, for any s greater than —1 and 
any a G and for any positive t, we have 

ct 

ii;)"T7. „, ,£/'+' .M|2 ^ r 



" JO I- h 



Proof. We shall only sketch the proof, as it is very close to the proof of Lemma 13.11 which was 
carried out above. The only difference is that the horizontal divergence of w_ does not vanish 
identically, but that will not change very much the estimates. We shall only write the proof 
in the case when a = and — 1 < s < 1, and leave the general case to the reader. Using 
Lemma 1.1 of [2j we have, for any in M, 

v^{t,-,y^)-Vhiif{t,-,y^)\uf{t,-,yi)) . < C||V,,i;(t, •, ys)!!^^ l|V,^2«^(t, •, ys)!!/^^ •, 2/3)11^- 

/ //^ h h 

Thus we get 

— Wml.^s + ||V,u;^(t)||^ , < \\\VHW%t)\\l,^^+C\\VHv!'mlo.M{t)\\l.Hs 

I IPi ,^ „. M„.,e-3/+ „, \\ j„, / l^h^ I. \\„.,e.hi 



-e' / {d:ip {t,-,y^)\vf-'\t,;yz))Hsdy:^- / iy P,{tr,y^)\w''\t,-,y^))^sdyz 
By integration by parts we have, thanks to the divergence free condition on uf, 
- / {d3PAt,-,y3)\lif'^{t,-,y3))H^dy3 = / (p, (t, •, y3)|932«^'^(t, •, y3))^s(iy3 



{p^{t,-,y3)\<iivhwf'''{t,-,y3))H-dys. 

By definition of the inner product of Hf^, we get 

{03PAt,-,y3)\uf'^{t,-,y3))Hsdy3 = / {VhpAt,-,y3)\uf'''{t,-,y3))H^dy3. 
I Jr 

Thus we have 

+ C\\\7hv\t)\\l^^,\\nf{t)\\l,^,-{l-e^) [ {V^pM, ■,y3M'\t, ■,y3))Hsdy3. 
Now we notice that 

2 

-{e'^di + Ah)p^ = div(/ • Whuf) = div^5j(i;V"), 

i=i 

which can be written in the simpler way 

-{£^dl + Ah)p^=diYhN'' 

with A^^ = ■ Vhuf'^ + d3{uf'^v^). It is easy to check that for any cr G M, 
(3.5) ||V;,pJ|^.^. <||iV^||^.^., 
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2a+4\ jCfh/c\\'2 



< ll/V'^ll^ 



We infer from p.Sp that 

l|V/.pJli2^. < \\i/''-Vhiif'''\\L2Hf: + \\Um''h'')\\L2H^ 

We claim that for all — 1 < s < 1, 

^h'^= / \C^hPAt,-,y3)\w'{t,-,y3))\fjsdy3 

(3.6) < \\\^hnf\\l2Hs+C\\wf\\l,^,\\VvX^Lli'^ + \\^^^^^^ 

Let us prove the claim. Suppose first that s = 0. Then a product law gives 

By interpolation, we get 

3 11 1 

The claim in the case s = follows from a convexity inequality, which gives 

Ih < \\\Vhvf\\l,^^,^ + C\\vf\\l,^^,^{l + \\v\\l^Ll)\\'^v\\l^Ll- 

In the case when < s < 1, then we use the product rule 
along with the fact that 

1h < ||V/iP 11^2 /f^-l II VfeU)^ 11^2 J^s. 

Finally in the case when — 1 < s < 0, we write 

< lll^'' • Vw^|L2fVs||lil^||r2rrs + ||W^ • U''||r2f>^ 

< l|V/||^ooi2||«;^|1^2^.||V,u;^|1^2^. + llli'^ll ||«;^|| +1 ||V,^^||^2^.. 

The claim (j3.6p follows by interpolation. 
Using that result we obtain that 

YM(t)\\l2H^ + l|V.2^^(i)||i2^. < l|V/(t)||i^^2||2£^(t)||^2^.(l + ll/(t)llioo^2) 
at ^vJ^h ^vJ^h ft ^yJ^h ft 
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and we conclude by a Gronwall lemma. Indeed we get that 



Ww'iml^H^ + [ \\yhwf{t')\\i,^^ dt' < \\wo\\i. 



X exp {c \\Vv\t')\\l^^,^{l + \\v\t')\\l^^,^ dt^ . 



But by the basic energy estimate (j3.ip . we have 

\\v!'{t)\\L^Li < II^^o|Il-L2. 
Moreover, Corollarv 13.11 implies that 

so Lemma 13.21 is proved in the case when q = 0. The case when a is positive is an easy 
adaptation of the proof of Lemma l3.lt it is left to the reader. □ 

Clearly Lemmas 13.11 and 13.21 allow to obtain Lemma 12.11 stated in the previous section. 

4. The estimate of the error term 

In this section we shall prove Lemma 12.21 stated above. We need to write down precisely the 
equation satisfied by the remainder term R^, and to check that the forcing terms appearing 
in the equation can be made small. 

Let us recall that 

= {{v!',Q) + e{'uf'^,e-^uf^^)){t,Xh,ex^) and 
Plppit,x) = {p^ + ep^){t,Xh,ex3). 

It is an easy computation to see that 

{dtVapp + Vapp ■ V<pp - Avlpp) {t, Xh, X3) = {dtV^ + v!" ■ W hV^ - Ahv!", 0)(t, Xh,eX3) 

+ e {^dtnf''' + v'^ ■ Vhnf''' - Aj^vf'"^ - e'dlnf'\ o) {t, xn^ex^) 
+ (0, dtnf'^ + / • Vhnf'^ - ^huf'^ - e^dluf'^^ (t, x/,, exg) + sF'it, Xh, exs) 
where 

F^t, Xh, ys) =^ {{euf ■ W'\wf ■ Vnf'^) + {vf ■ V^;^ 0) + e{dlv\ 0)) (t, x^, ys). 
In order to simplify the notation, let us write F"^ = F'^'^ + F^'"^ with 

F^'i # {ew' -Vw'^^^w' ■Vw!^^) + {'uf -Vv^^Q) and 

F-'2 def ^(02^^^ 0), 

Recalling the equations satisfied by and nf, we infer that 

{dtvlpp + vlpp ■ Vvlpp - Avlpp) {t, Xh, xs) = -Vplpp + eCit, xt, exs) 
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with G^{t,Xh,yz) =^ (^F^ + (0,53Pg)) {t,Xh,yz) and F^{t,Xh,x^) eG^{t,Xh,ex^). Denot- 
ing = — p%pp., we infer that 

dtR" + ■ VR" + vlpp ■ VR^ + R" ■ Vvl^p - AR" = -Vq^ + F^. 

So Lemma \T2\ will be established as soon as we prove that „, . ■ i , < C„„ ^„e3. 

The forcing term consists in three different types of terms: a pressure term involv- 
ing p^, a linear term e'^d^v^{t,Xh,£X3), and finally a number of nonlinear terms, defined 

as eF^'^{t,Xh,£x^) above. Each of these contributions will be dealt with separately. Let us 
start by the pressure term. 

Lemma 4.1. The following estimate holds: 



Proof. We define Pq (t, xs) {d^p ){t,Xh,£Xs). Sobolev embeddings enable us to write 



-0' 

^ 1 1 T~tF 1 1 

2 , 



' ° "L2(R+;ij~i{R»)) ^ " "l2(R+;L2(r3)) 



Recalling that 

2 

we have by Sobolev embeddings, 

2 



< 



1 1 1^ 1 1 L°° (R+ ;L3 (R3 ) ) 1 1 (^3:?i 1 1 L2 . (ig3 ) ) 



5, ll^ll j_ 1 . .J.. Iio'st'll , j_ 1 

~ II — 1 1 r 9 ^Tre+ . rr?T /m)3\\ 1 1 — 1 1 r nr-. /m3+ . rr^r / 



'L2(R+;ii'^(R3)) " '^-"L°°(R+;_H'^(R3))' 

SO we can conclude by Lemma |3. 11 This proves Lemma l4.ll □ 

Now let us consider the linear term e^c?|w^(t, x/i, exs). The statement is the following. 
Lemma 4.2. The following estimate holds: 

Proof. We have 

e'^\\(d'^v^)(t,xh,ex3)\\ , , ■ i < e d3(d3V^(t,Xh,ex3)) i 



< emv''){t,Xh,X3)\\^ 1 
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A computation in Fourier variables shows that, for any function a on M'^, we have 



< 



1 



1 



By interpolation, we deduce that 

|2 



6 



d^3 



\Ch\dCh—+ , 



6\ 2 161 d^3 



3^||V,,a|| 

Applying this inequality with a = d^v, we get 

1 



l2(ir3) + l!a|lL2(]g3)||93a||^2(]R3). 



L2(IR+.i2 



\dhit)f 



1 

2 

L2{R+;L2 
1 



L2(R+;L2 



by Lemma l3.ll This proves Lemma l4.2[ 



'£,1 



□ 



Now let us turn to the nonlinear terms composing F'^, which we denoted above eF' 
Lemma 4.3. The following estimate holds: 

e||F^'i(t,x,,ex3)||^,^^^^^_i^^3^^ < Cy,,y,,el. 

Proof. We recall that 

Notice that for all functions a and b and any 1 < j < 3, 

II°^^-^IIl2(M+;H-^(R3)) ^ ll"^^-^llL2(R+;Li(R3)) 

^ ll"llL°°(R+;L3(R3))ll^i^llL2(K+;L3(R3)). 

Defining c^(t, x/i, 3:3) = {adjb){t,Xh,£X3) this implies that 

_ 2 

We can apply that inequality to a and b equal to t; or y;^, due to the results proved in Section [21 
and the lemma follows. □ 
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